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Effective quantum field theoretical continuum models for graphene are investigated. The models 
include a complex scalar field and a vector gauge field. Different gauge theories are considered and 
their gap patterns for the scalar, vector, and fermion excitations are investigated. Different gauge 
groups lead to different relations between the gaps, which can be used to experimentally distinguish 
the gauge theories. In this class of models the fermionic gap is a dynamic quantity. The finite- 
energy vortex solutions of the gauge models have the fiux of the magnetic field quantized, making 
the Bohm-Aharonov effect active even when external electromagnetic fields are absent. The fiux 
comes proportional to the scalar field angular momentum quantum number. The zero modes of the 
Dirac equation show that the gauge models considered here are compatible with fractionalization. 
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I. INTRODUCTION AND MOTIVATION 

In a monolayer of graphenet"— , the singlc-particlc dis- 
persion relation near the so-called K and K' Dirac points 
is linear in |fc| - see, for examplcf^ii. Formally, it is the 
dispersion relation of a massless relativistic fermion. Fur- 
thermore, the description of the low energy electronic 
excitations can be accommodated in a Dirac-type equa- 
tion. Indeed, starting from a tight-binding Hamiltonian, 
with a hopping parameter independent of the lattice site, 
one can compute exactly the dispersion relation, expand 
around the two inequivalent Dirac points to rewrite the 
dynamical equations of motion for electrons and holes as 
a Dirac equation in two dimensions with the four com- 
ponent spinor given by 

The indices a and b refer to the two triangular sublat- 
tices, while the -I- and — indices to the two Dirac K 
and K' points. For a perfect graphene crystal structure 
the fermions behave as massless relativistic particles^, 
which translates into the well known ballistic behavior 
of the electrons2r2ij and there is no gap between va- 
lence and conduction bands. However, if the two dimen- 
sional honeycomb array of carbon atoms is distorted due 
to the presence of impurities or to the distortion of the 
crystal structure, for example, the fermions acquire an 
effective gap given by half of the mass gap. Besides 
fermion mass generation, quantum Hall effedj^iiSrii, 
fractionalizationi^"— and Berry phasesi^ have been ob- 
served in two dimensional graphenelike structures. 

Fractionalization in one dimensional models was inves- 
tigated more than three decades ago^"— within poly- 
acetylene. A similar phenomena like fractional quantum 
hall effect due to quasiparticle fractional charge and/or 



fractional statistics can take place in two dimensional 
systemsi^ii^iii. 

A dynamical theory for two dimensional graphene 
should describe, of course, its phenomenology and should 
be able to accommodate for the possibility of gap gen- 
eration and fractionalization. Ini^ the authors presented 
a mechanism for electron fractionalization in graphene- 
like systems keeping time reversal symmetry. Invoking a 
Kckulc texture, a complex order parameter Aq was in- 
troduced. Aq couples the two Dirac points and changes 

the electron dispersion relation to e(fc) = \Jk?-\- | Aop. 

Assuming a vortex-like profile A(r) = A(r)e™^, where 
n is an integer, with A(r) vanishing as r'"' for small r 
and approaching Aq at large r, fractionalization was as- 
sociated with the presence of a zero mode of the Dirac 
kernel. Although fractionalization was connected with 
the vortex shape of A(r), vo^ the authors do not specify 
the dynamics for the complex vortex profile. 

Ir*i^ a dynamical content to the vortices was intro- 
duced through a chiral gauge theory which is compati- 
ble with electron fractionalization, extending the work of 
Hou et al^. In the language o&^, vortices are associated 
with a complex scalar field ip which couples linearly to the 
fermions. Although, a dynamical equation was written, 
see equations (13) and (14) ini^, the potential V{(p*ip) 
was left unspecified. From the point of view of Quantum 
Field Theory (QFT), there is no reason to exclude other 
types of symmetries and couplings, not present in the 
model discussed by Jackiw and Pi, without destroying 
fractionalization. 

Indeed, field theoretical models have been applied 
to describe nanotubes and graphene physics with some 
success in reproducing their quantum properties (see 
e.gi^i^). In this work, we elaborate on derivative free 
fermion-boson and boson self interactions allowed by 
QFT principles and discuss possible gauge interactions. 
The models considered here are a generalization of the 
results ofi^ and, besides the fermionic field, they con- 
sider a complex scalar field ip and a single gauge vector 
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Af^ field. Moreover, possible ways to distinguish between 
the different gauge symmetries are discussed. 

Our interpretation for the complex scalar field and 
gauge field being that ip and resume all the dynamics 
of the self-interaction of the carbon background and the 
mean fermionic self interaction. 

The usage of scalar and vector potentials to describe 
some of the graphene properties is not new. Indeed, 
scalar and vector potentials, including gauge fields, have 
been used in the literature to describe disorder phenom- 
ena, including distortions of the lattice honeycomb, struc- 
tural defects, point defects and self-doping effects asso- 
ciated with the breaking of electron-hole symmetry near 
the Dirac points among other properties. A detailed dis- 
cussion can be found inA^S^ and references there in. 

Graphene is an electrical neutral system. On the other 
hand graphene has charge carriers. Therefore, it seems 
natural to associate a charged field with the carbon back- 
ground. Furthermore, if resumes the carbon back- 
ground self-interactions it should be able to accommo- 
date for the propagation of phonons in the carbon lat- 
tice. Phonons feel the density of states of the fermionic 
degrees of freedom and one expects if to couple to the 
density of electron/holes, i.e. to iptp ~ — (V-'+)^'0- ~ 

V- - {i'-Y^i - (^Z--) V+ • Throughout this pa- 
per, we will use the chiral representation for the Dirac 
matrices, where 

and 

-^(j"-0^ (3) 

stand for the Pauli matrices. Besides the coupling to 
the density of electron/holes, QFT allows also for a pseu- 
doscalar like interaction, described by a coupling of tp to 

The scalar and pseudoscalar interactions couple the two 
triangular sublattices and the two Dirac points K and K' 
in different ways. The models discussed in the present 
work explore contributions coming from both type of in- 
teractions (i.e., ipip and ipj^ip)- 

The potential energy for the complex scalar field ip can 
accommodate a non-vanishing vacuum expectation value 
{(fi). If {(fi) 0, then the model generates a fermion mass 
via spontaneous symmetry breaking. On the other hand, 
if ((/?) = 0, the electrons in graphene remain gapless. 
Therefore, we identify pure graphene with the vacuum 
state where {(p) = 0, with all other graphene distorted 
and/or doped states begin described by different vacuum 
and for these (cp) ^ 0. 

In what concerns the bosonic fields, the model can ac- 
commodate mass gaps both for the scalar and vector ex- 
citations. In general, for the gauge theories considered 
here, a fermion mass gap implies also a vector mass gap. 



This comes directly from the Higgs mechanism for mass 
generation. The gap for the scalar excitations is linked 
with the details of the potential V{(p'^p) and is not di- 
rectly coupled with the fermion and vector gaps. Indeed, 
we found that the scalar gap can vanish independently 
of the fermion and vector gaps. 

In this paper we also discuss a number of differ- 
ent gauge models which, in principle could be suit- 
able to describe graphene properties. The relation be- 
tween the spectrum of the scalar and vector excitations 
with the fermionic spectra depends on which symmetry 
is gauged. Further, the different connections between 
fermion, scalar and vector gaps opens for the possibility 
to check experimentally which of the gauge symmetries, if 
any, is realized in graphene. Changing the fermion mass 
gap, for example modifying the concentration of impu- 
rities and or the distortion of the lattice, and looking 
on how the scalar/vector mass excitations adjust them- 
selves, one can distinguish between the different gauge 
models. Besides the pattern of the mass gaps, in general, 
the models allow also for vortex like solutions and are 
compatible with fractional statistics. 

These gauge models have finite energy vortex solu- 
tions. For one example, we show that the vortex solu- 
tion implies in the quantization of the "magnetic field" 
fiux. In this case, the abelian gauge field is connected 
with the angular momenta of ip along an axis perpendic- 
ular to the graphene sheet and, in this sense, the fiux 
of the "magnetic field" is a measure of the angular mo- 
menta of ip. This particular solution can be interpreted 
as consequence of topological defects in the graphene 
structure and, in principle, phenomena like the Bohm- 
Aharonov^"— effect can occur even when external elec- 
tric and magnetic fields are absent. 

We show that the vortex solutions of non-chiral gauge 
models presented here, have normalizable zero modes of 
the Dirac equation. The presence of the normalizable 
zero modes implies fractionalization for graphene and 
the quantum Hall effect in graphene sheets could become 
possible even without chiral gauge symmetry and with- 
out external electric and magnetic fields. The observa- 
tion of the quantum Hall effect in two dimensional ma- 
terials without external electromagnetic fields was also 
discussed, within the framework of tight binding mod- 
els, in the work of Haldane^S and Hill et al^. Accord- 
ing to the later work, the observation of the quantum 
Hall effect without an electromagnetic field requires the 
breaking of the sublattice symmetry, where the two sub- 
lattices a and b are interchanged, and the opening of a 
mass gap at one of the Dirac points, let us say K, while 
the other Dirac point K' remains gapless. In the class 
of models discussed here, the mass gap is open, or not, 
simultaneously at K and K' . Then, according to^ quan- 
tum Hall effect without an external magnetic field is not 
measurable as the contributions from K and K' to the 
Hall conductivity cancel exactly. The models considered 
in the present work, although reproduce the tight bind- 
ing model in the appropriate limit (see, for example, the 
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paper—) they give dynamics to all the fields that repre- 
sent the carbon graphene background and the excitations 
tf and the gauge field. Recall that is charged and can 
give rise to an electric current. In this sense, the models 
go beyond the tight binding model, opening the possibil- 
ity of having a dynamical situation where the conditions 
explored ioSi do no apply, and, perhaps, may allow the 
measurement of the quantum Hall effect in graphene. We 
plan to address this question in a future publication. 

We would like to call the reader attention that, within 
the class of gauge models discussed in the present work, 
fractionalization is allowed without breaking any of the 
usual discrete symmetries like, for example, time reversal. 
We do not compute the rich set of zero modes of the 
Dirac equation but those ones obtained here are, again, 
connected with the angular momentum of the complex 
scalar field tp. 

This paper is organized as follows. In section |ll] the 
effective QFT for graphene are discussed, including its 
global symmetries. In section IIIIl the different global 
symmetries are gauged and we discuss how this change 
the spectra of the scalar, vector and fermion excitations. 
Furthermore, combining the information on the different 
types of gaps, we are able to suggest an experimental 
test to disentangle which of the gauge symmetries ap- 
ply to graphene. In section IIVI the equations of motion 
are derived and the vortex solutions for the gauge mod- 
els are discussed. The short distance and long distance 
properties of the vortex are computed explicitly. The 
gauge models predict the fiux quantization of the "mag- 
netic field" associated with the gauge field. Further, the 
flux quantization is connected with the angular momenta 
of <p. In section |V] the zero mode solutions of the Dirac 
equation for a vortex configuration are investigated. Fi- 
nally, in section IVTl we resume and conclude. 

II. THE EFFECTIVE MODEL 

Let us assume that the charge carries, i.e. electrons 
and holes, are relativistic fermions described by a four 
component spinor ij). The Lagrangian density describing 
the interaction between fermions ip and (p can be written 
as 

-P{^)^^ ~ P5M?^75V^- V{p^ip), (4) 

where the polynomial P{^p) and P5((/?) define the type 
of interaction between fermions and the carbon crystal 



structure and V{ip^'p) the self interactions of the back- 
ground structure. 

The reader should note the linear combination of P tpip 
and P5 '075'0 couplings. Such a freedom would allows to 
set different couplings to each of the possible fermion chi- 
ralities and, in this way, build a chiral theory. Moreover, 
other Dirac 7-matrices are allowed. However, to keep it 
as simple as possible and to avoid derivative couplings, 
in the following we will consider only scalar and pseu- 
doscalar types of interactions. 

In a system of units where the action is dimensionlcss, 
space and time have dimensions of an inverse mass, and 
for two spatial dimensions and one temporal dimension, 
"0 has dimension of mass, [ip] ^ M, and ip of square root 
of mass, [ip] ~ M^/^. Requiring that the theory described 
by C is perturbatively renormalizable, then for polyno- 
mial interactions, naive power counting forbids coupling 
constants [g] M" with a < 0. This fixes unambigu- 
ously the interaction terms to 

P{f) = 91 if + f'') + 92 </?V , (5) 

P^ip) = hi (p - p^) + ih2ip^ip, (6) 

and 

Vip^p) = ^,^pW) + y(^V)' (7) 

up to a constant Vq. 

If in P{p) and P^iif) one takes hi = —gi and 52 = 
/i2 = one recovers the Jackiw-Pi theory with their p"" = 
2 Re(i^) and p^ =2 Im((^) ~ see equation (8) ini^. In this 
sense C generalizes the results ofi^. 

Let us discuss now the global symmetries of the model 
described by the Lagrangian density (|4]) . 



A. Global (7a (1) Symmetry 

One of the motivations ofi^ was to build a chiral gauge 
theory. So let us consider the same type of chiral trans- 
formation, i.e. 

^ — y e^^i^ <p ^ e*" . (8) 

To first order in w and 77, the corresponding variation of 
the Lagrangian density reads 



I 

|2ia; [hi (p - p^) + z/i2 (y'V)] + iV9i (</' - <y3^) } -^^A + {siw [gi {p + p^) + 52 (vV)] + (</' + } V^^sV' ■ (9) 

I 

Requiring invariance of £ under the transformation it di = i^i and ?/ = ±2w, with the minus sign recovering 
follows that g2 ~ h2 — 0, as in the Jackiw-Pi theory, and 
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the original Jackiw-Pi theory. Note that, from the point 
of view of the ([5]), invariance of the theory means that 
the chiral charge associated with Lp is, up to a sign, twice 
the fermionic charge. 

The set of transformations (|8]) witli ?/ = ±2w form a 
group which will be called from now on J7a(1)- Recall 
that the Lagrangian density is invariant under ?7a(1) if 
and only if 32 = ^2 = and g\ = zLhi. 

B. Global C/(l) Symmetries 

Besides the chiral transformation just discussed, the 
Lagrangian density (j4|) has further non-chiral U{1) global 
symmetries. The set of transformations 

ip — ^e'^V (10) 
defines the Uf{l) global symmetry of C and the set 

(^-^e'^V', (11) 

defines the Ub{l) global symmetry of £ if gi = /ii = 0. 
Further, ii gi ^ hi = in ([S]) and ([B]), then the set of 
transformations 

iP — >e"^il;, if — ^e'^'ip, (12) 

where uj and 77 are independent parameters defines an- 
other global symmetry of £, called below J7/(l) (X) [/f,(l). 
Gauging this symmetry requires the introduction of two 
gauge fields unless one imposes an additional discrete 
symmetry with respect to the interchange of the gauge 
fields. Note that the discrete symmetry gives no con- 
straint on the coupling constants for the fermionic, g, 
and bosonic, g^p^ fields. Indeed, the fermionic covariant 
derivative reads = d^^ + igA^, while the bosonic co- 
variant derivative is = + ig^A^^, where A^^ is the 
gauge field. For the sake of simplicity, i.e. to avoid con- 
sidering more than one gauge field, we will analyze only 
the Uf{l) (E) Ub{l) symmetry supplement with the dis- 
crete symmetry. Anyway, we will keep using the name 
Uf{l) ® Ub{l) for the symmetry. 

C. On The Various Global Symmetries 

The various global symmetries are distinguished by the 
nature of the (p — tp interaction and by the number of in- 
dependent coupling constants associated with the gauge 
field required to define the model. 

The J7yi(l) symmetry is not compatible with the in- 
teractions {f^f) V'V' s-nd {(p''(p) ipjsip s-iid only linear 
in ip terms are allowed in the interaction with fcrmions. 
Further, the model defines a unique coupling constant. 

The Uf{l) allows linear and quadratic ip couplings to 
the fermionic field and requires a unique gauge coupling 
constant. 

The local symmetries 14(1) and Uf{l)(E)Ub{l) are com- 
patible only with a quadratic (p coupling to the fermion 



Symmetry 


Constraints 




.92 = h2 


= 0, gi = ±hi 


Uf{l) 


none 






51 = hi 


= 


[7/(1)® [/,(!) 


gi = hi 


= 



TABLE I: Global symmetries of the Lagrangian density Q. 
Recall that for [7/(l)(g)C/i,(l) we impose an additional discrete 
symmetry - see text for discussions. 

fields. If Ub{l) requires a unique gauge coupling constant, 
the gauge model with Uf{l)<^Ub{l) as a symmetry group 
includes two independent gauge couplings. 

Tabic H] resumes the global symmetries of Lagrangian 
density (U) and the corresponding constraints on the ip—tp 
coupling constants. 



III. GAUGE MODELS AND MASS GAP 

The various U{1) symmetries of C, see table U can 
be made local. Different symmetries will lead to dif- 
ferent gauge theories for graphene, after replacing the 
derivatives by covariant derivatives and adding the cor- 
responding kinetic term for the gauge field. Naturally, 
the different symmetries will introduce different dynam- 
ics which can be seen, for example, at the level of the 
theory spectra, i.e. at the various mass gaps. 

Recall that we are excluding derivative type couplings. 
In what concerns the gauge field, not including derivative 
couplings means that a priori we are excluding a Chern- 
Simons ternt^"— 

e"'5^^„ {dM (13) 

in C. This type of interaction is allowed by gauge in- 
variance and, in 2-1-1 dimensions, is not excluded by the 
renormalizability requirement. 



A. Scalar Mass Gap 

We start our discussion looking at the scalar excita- 
tions in graphene, i.e. looking at the mass spectra for 
the complex scalar field tp. 

The self interactions of (p are described by the potential 
energy ¥{1^^^) ~ see equation Depending on the val- 
ues for /i^ . A3 and Ag , V can have cither one, two or three 
minima. The discussion of the V{ip^(p) cxtrema is rela- 
tively straightforward and will not be reproduced here. 
The relation between potential parameters and number 
of extrema is summarized in table |lll 

The mass gap for the scalar excitations, i.e. the mass 
associated with the complex scalar field, can be computed 
from I?]) writing tp ~ v + where v = {tp) is the vacuum 
expectation value of ip assumed to be real. If ip can be 
rotated in such a way that it becomes a real field, then 
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Item 



2 

M 


A4 




jr 



^ extrema 


< 


> 


< 





1 maximum 


> 


< 


> 





1 minimum 


> 


> 


> 


> 


1 minimum 


> 


< 


> 


> 


5 extrema 


< 


any 


> 


> 


3 extrema 



TABLE II: The number of extrema of V{^) as a function of 
the potential parameters. Our definition for 5 is S = \\ — 

it follows that the quadratic term in is given by 

i Ml $2 = (A4 + Ae) $2 (14) 

and one can define the mass gap for scalar excitations as 
A$ = = 2 \v\ v/A4 + 2v2A6 . (15) 

The scalar mass gap is, then, independent of the gauged 
symmetry. 

A non-vanishing Af$ requires either a non-vanishing 
expectation value for {ip) or a {i^) = and a, fj? > 0. Fur- 
thermore, a no n- vanishing scalar gap requires also that 
A4 + 2u2A6 > 0, if {if) ^ 0. 

From equations ([Til) and ([T51) it follows that the model 
is compatible with a non- vanishing vacuum expectation 
value, i.e. a {p) ^ 0, in combination with a vanishing 
scalar gap if and only if A4 = —2v^Xq. In this case, the 
theory predicts a fermionic mass gap proportional to {tp), 
with no gap formation for the scalar excitations, i.e the 
dispersion relation for the scalar excitations is linear in 

1^1- 

B. Fermionic Mass Gap 

For the fermions fields, if tp and/or 'p'^p acquire a non- 
vanishing vacuum expectation value, then C acquires a 
mass term and a chiral mass term - see equations Q, 
© and ©. 

Let us assume that {ip) = u 7^ 0, with v being a real 
number. From the point of view of the fermions them- 
selves, the interaction with the carbon structure shows 
up as 

fhipip + j /i2 75 "0 (16) 

where fh = 2 giv + g2V^- There is no reason a priori to 
require the positivity of m or /i2- Indeed, solving the free 
Dirac equation, with a mass term given by (|16p , gives the 
following dispersion relation 

e(p 2) = ^p^+m^ + hlv^, (17) 

i.e. the effective fermion mass is given by -y/ + h'^v'^ 
and the corresponding mass gap between valence and 



conducting bands is 2 ^Jrh? + h^v^ . The above reason- 
ing is valid even when {^p) — and {ip^ip) 7^ 0. In this 
case TO — g2{(p^p) and h2{(p^(p) replaces h2V^. 

The coupling of the complex scalar field tp to the 
fermion degrees of freedom is given by 

(51 {^ + ^^) + 52 (^V) iptp 

+ (^hi{(p-p^) + i/i2V5V) V'TsV', (18) 

which has exactly the same structure as the mass term 
given by equation ((T6)) . Therefore, the status of the field 
ip can be translated into a dynamical fermion mass, i.e. 
a dynamical mass gap, which is both time and spatial 
dependent. The model accommodates graphene states 
where for certain space-time regions the system is gap- 
less, i.e (p 0, and for others where (p = and there is no 
gap. We are currently exploring the implications of this 
dynamical gap to graphene properties and will report the 
results elsewhere. 

In graphene the fermionic mass gap A is a function of 
gauged symmetry. It follows that for 

Ua{1) A^A\giv\, (19) 

Uf{l) A = 2^?h? + hlv^, (20) 

C/fc(l), (7/(1)® C/fc(l) A = 2^gl + hl (21) 

where A is twice the fermion mass. 



C. Vector Mass Gap 

It remains to discuss the mass gap for the vector ex- 
citations in graphene. The mass term for is gener- 
ated by the scalar kinetic part of £, i.e. by {p^l^p>)^ D^ip. 
Therefore, unless the gauge transformation changes ip, 
the gauge field remains massless. It follows that for 

Ua{\) AA^^\giv\, (22) 

f//(l) Aa=0, (23) 

Vb(X),Vf{\)®Ui,(X) AA^^f2\g^v\ (24) 

where A^ is the vector mass, i.e. the vector mass gap. 

D. Gauge Symmetries and Gap Relations 

The mass, i.e. the gaps, for each of the fields in the 
model are generated via Higgs mechanism. Besides the 
mass, the Higgs mechanism also provides a relation, de- 
pendent on the symmetry group, between the different 
gaps - see Table IIIII for a summary of the results dis- 
cussed in the previous sections. 

A U AiX) chiral gauge theory implies a linear relation 
between the vector and fermion mass gaps A^ = A/ -^/S, 
while C/h(l) or J7/(l) ® UhiX) relate the two mass gaps 



6 



Symmetry 


A 


Aa 


Ua{1) 


4 \9iv\ 


V2 \gi v\ 


Uf{l) 











V2 \g^v\ 


Uf{l)(SU,{l) 




V2\g^v\ 



TABLE III: Mass gaps as a function of the gauge symmetry - 
see text for discussions. The scalar mass gap is independent 
of gauge group and for a non-vanishing vacuum expectation 
value if given by A$ — 2\v\ \'\4 + 2v^X^. 




A 



FIG. 1: Aa and A# as function of the fermionic mass gap A. 
The curves where computed setting all the coupling constants 
and the potential parameters to one and using arbitrary units. 



by a quadratic equation A\ = A g^j \J g\ + For 
the gauge theory associated with the J7/(l) symmetry, 
there is a mass gap for the fermionic and scalar degrees 
of freedom, but no gap for the vector excitations. 

A non- vanishing fermionic mass gap requires a (ip) ^ 
or {(p'^if) 0, which by itself implies a scalar mass gap in 
graphene. The connection between the scalar mass gap 
Aip and the remaining gaps is slightly more complicated 
than the relation between A and A^. 

The scalar, fermion and vector mass gaps arc func- 
tions of the coupling constants and of the (p vacuum ex- 
pectation value. The connection between the mass gaps 
and the fundamental parameters of the theory depend 
on which global symmetry is gauged, as summarized in 
table Hill Moreover, all the mass gaps can be written in 
terms of the fermionic gap A. Therefore, if one is able 
to build graphene with different mass gaps, for example 
changing its doping and/or distortion, one can test which 
of the symmetries discussed before applies to graphene, 
if any, simply looking at how Aa and A* change with 
A. As an illustration, in figure [1] we show A^ and A$ 
as a function of A when all the coupling constants are 
set to unit. The figure uses arbitrary units. Clearly, the 
functional behaviour distinguish between a chiral gauge 
theory and a Ub{l) or a C//(l) (8) gauge symmetry. 

A Aa = for all values of A clearly points towards a 



(7/(1) gauge theory. 

Of the gauge theories discussed here, the mass gap re- 
lations do not distinguish between the two symmetries 
Ub{l) and Uf{l) (E) Ub{l) as they provide similar types 
of — Af^ interactions. However, if C/f,(l) does not cou- 
ples directly the fermions to the vector excitations, the 
Uf{l) (g) Ub{l) gauge symmetry has such a direct cou- 
pling between electrons and gauge field. Of course, if 
the strength of this coupling is extremely small, then the 
two theories will give exactly the same predictions for 
graphene. However, if the coupling tp^^^ip is sizable 
enough, then the fermionic interaction with the vector ex- 
citations will give relevant contributions to the dynamics 
of graphene and the two symmetries will provide different 
physics for vector and fermion excitations. 

In the models discussed above, the gap in graphene is 
generated via spontaneous symmetry breaking. This is 
not the only way of having massive particles in a theory. 
For example, as discussed iri^, a fermionic gap can also 
be generated via dynamical symmetry breaking. 



IV. EQUATIONS OF MOTION, VORTICES AND 
FLUX QUANTIZATION 

The equations of motion associated with the various 
fields are derived from C in the usual way. For fermions, 
they are given by 

tj^D^^ - P{^)^ - PsMts^ = 0, (25) 

where = 9^ + igA^i is the covariant derivative and A^ 
is the gauge field. The corresponding equation for Lp is 

Df^D^ip ^ - (51 + 92 v) ^V' + (^1 - ih2 v) V' 75V' 

- H^ip - A4 (((sV) "i^ - Ae (((sV)^ </' (26) 

with the covariant derivative given by = igip^p,- 
The gauge field equation of motion reads 

d^F-^'^ ^ ~g^j''yj-tg^pHD''p) + ig^{D''p,)^ p. (27) 

Equations (j^S]) , ((25)) and ([?7)) are the equations of mo- 
tion derived from C taking into account all the possible 
coupling constants. In order to study each of the gauge 
theories considered previously, one has to take into ac- 
count the corresponding constraints associated with the 
gauge group and summarized in table HI 



A. Vortex Solutions 

Let us now discuss vortex like solutions for the bosonic 
sector of the theory, disregarding the coupling to the 
fermionic degrees of freedom. This provides a consistent 
solution for the field equations for the Dirac zero modes 
discussed in the next section. 
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In this section wc will consider static solutions for equa- 
tions (HH) and (HZl) with 



and 



where 



0, 



A' 



a'{r) 



da{r) 
dr 



bir). 



(28) 
(29) 

(30) 



The equation of motion for (f then becomes 

= 

^ fi^ipo + Xiifil + Xeifil (31) 
and the gauge field equation simplifies into 



1 d 


' dipo' 


ni 




r dr 


dr 


V r 





1 d 

r dr 



db 

dr 



2g^ C-^+g^b^^l. (32) 



Note that only for C/^(l), Ub{l) and t//(l) (g) C/f,(l) gauge 
theories, i.e. when g^ ^ 0, the scalar and gauge fields are 
coupled. In this section we will consider only the Ub{l) 
and [//(I) Ub{l) symmetries. 

Let us first discuss the solutions of equations (|3T|) and 
([5^ at small distances. For r <^1, the gauge field equa- 
tion becomes 



l_d 

r dr 



db 
dr 



= 0. 



(33) 



provided that ipo{r) is regular at the origin, whose solu- 
tion is 



b{r) = fei r H 

r 



(34) 



where &i and 6_i are constant of integration. If &_i = 0, 
equation (PT]) 



l_d_ 
r dr 



d(£o_ 

dr 



and its power law solution reads 
Mr) = 



(35) 



(36) 



On the other hand if bi = 0, in ((35)) should be re- 
placed by (n-f g^6_i)^ and the corresponding scalar field 
solution at small r is 



ipo{r) = r 



(37) 



It follows that ipoif) is always regular at the origin. 

At large distances, for finite energy solutions, ip ap- 
proaches a constant value, a minimum of V{iPq), and the 
l.h.s of equation (j3ip vanishes. If (po ^ 0, then 



bir) 



n 

g^r 



(38) 



and b{r) is also a solution of equation ([32]) . Note that this 
solution can be extended to full range of r values, with 
the exception of the origin, where a delta function sets 
in coming from the laplacian of b. On the other hand, if 
LPQ = 0, i.e. for pure graphene, one still has a solution for 
b as in ([34|). However, the requirement of finite energy 
demands bi — Q. 

The vortex solutions include, as a particular case, the 
type of configurations considered ir>i^, where b{r) is reg- 
ular and '^air) ^ r'"', for r ^ 1, and, for large r, (pQ 
becomes constant and b{r) ^ 1/r. A class of vortices 
whose short distance behavior is given by I\i7\ was found. 
Furthermore, the zero modes of the Dirac equation com- 
puted in the next section require a vortex solution with 
LpQ constant and non- vanishing and b{r) = oi'{r) given 
by (p8)) over all space, with the exception of the origin as 
discussed previously. From the point of view of the gauge 
models, the singular behavior at r = does not raise any 
conceptual problems. Indeed, we are using a continuous 
model to describe graphene and the dimensions of the 
unit cell provide a natural short-distance cut-off below 
which the model is no longer valid or, at best, should be 
corrected to take into account the crystal structure of the 
carbons and lattice defects. 

If one takes the usual definition for the "magnetic 
field", S = V X A, it follows that B vanishes at large 
r. Close to the origin, B approaches a constant for type 
(l36l) solution and vanishes for ([37]) configurations. For 
both type of configurations, the vortex energy 



d^x 



Dip 



Viip^ip) 



(39) 



is finite. 



Despite the vanishing of B at large distance, the "mag- 
netic flux" of the vortex configuration over a suffi- 
ciently large radius closed surface is quantized. Indeed, 
for a spherical surface centered at the origin. 



B-dS 



A-dl 



2tt r b{r) 



2n7r 

gv 



(40) 



where n = 0,±1,... is the component of the angular 
momentum on an axis perpendicular to graphene plan 
associated with the complex scalar field (p. Flux quanti- 
zation opens the possibility of having Bohm-Aharonov 
type effects in graphene without external electromag- 
netic fields, where electrons are scattered by a vector 
potential, which can be associated with a "topological 
defect" on the carbon structure, and acquire an extra 
phase. We call the reader attention that Bohm-Aharonov 
phases have been observed in suspended graphene, in as- 
sociation with mesoscopic deformations, where the mea- 
sured charge carriers mobility is substantial larger than 
in graphene on a substrate - seo2^ and references therein. 
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FRACTIONALIZATION, DIRAC EQUATION 
AND ZERO MODES 



As discussed at the beginning of the present work, elec- 
tron fractionahzation is related to the normalizable zero 
modes of the Dirac kernel. The presence of these zero 
modes opens the possibility of observation of fractional 
quantum Hall effect in graphene. A nice discussion con- 
necting the Dirac equation zero modes with electron frac- 
tionahzation can be found, for example, ii>i^. 

For the C/4(l) gauge theory, the proof of the pres- 
ence of normalizable zero modes can be found ini^. The 
Uf{l) gauge theory includes, as a particular case, the 
Dirac equation discussed by Hon, Chamon and Mudry 
ini^. Therefore, for Uf{l) gauge theory fractionahza- 
tion is possible. It remains to discuss the Ub{l) and 
Uf{l) CS) Ub{l) gauge theories. In the following, it will 
be assumed that the bosonic sector is in a static vortex 
configuration with A*^ — 0, 



= e^^dja{r) and (p{f) = (poir) 



(41) 



For a general gauge field, the Dirac equation associated 
with the t/f,(l) and Uf{l)®Ub{\) gauge theories is given 
by 

I - ia • (V - igA) 

+ .92 (<pV) + (l^V) ;^75|V' = E^j}. 

(42) 

Let us define the following function 



l\{r) = zipl{r) e 



92+^2^ 



where 



tana 



92 



(44) 



With the above definitions and for the Dirac spinor 



the Dirac equation becomes 



(45) 



idr - — - iga' ) + A(f) ^'L ^ E 



-idr + — + iga' ) + A* (f) = E '^t , 



-idr + y+ *.9a') + Mr) = E-^l, 

-idr-y-iga'^'ift+A*{f)^l = 

(46) 

where a' means the derivative with respect to r of func- 
tion a{r). These equations are invariant under the in- 
terchange of the two sublattices a -O- & provided that 
6^—9 and a' — > —a'. This symmetry generalizes the 
sublattice symmetry of the Dirac equation already dis- 
cussed ir^i^. We proceed assuming that = 0. Note 
that, for zero modes, given a solution of the Dirac equa- 
tion in sublattice a, the generalized sublattice symmetry 
generates another zero mode but leaving in sublattice b, 
or vice- versa. If one writes 



the zero mode equations become 

+ (7+50') (l)+ + z^l{r)4>. 



(47) 



(48) 
(49) 



(43) if the following relations 



fc = TO - 1 and /3+ = - a (50) 



are satisfied. From equation pS)) one can write 



Ao(r) 



— + .g a ) (p+ 



(51) 



where Ao(r) = z '~p\(r). Replacing this expression for 
in equation (|49p one arrives at the following second order 
differential equation 



Ao 



m 

ga ' - \ — 



ga 



k , A'„ 
-+ga' + -^ 
r Ao 



--A2 



0. 



(52) 



The computation of a solution of equation (|52[) requires the knowledge of (poir) and a'{r). Let us look for config- 



9 



urations where (po(r) is a non vanishing constant, i.e. a 
minimum of y((^g). Then, Ag = and Aq = zcpg is also 
a non vanishing constant. The equation of motion of the 
scalar field ([5T|) gives 



6(r) = a'{r) = 



(53) 



This particular gauge configuration solves the equation of 
motion for the gauge field ([32l) , except at the origin where 
a Dirac delta function sets is due to the laplacian opera- 
tor. Our vortex solution requires a short-distance cut-off, 
which is provided by the dimensions of the graphene unit 
cell or the length scale associated with a defect. Indeed, 
for distances smaller than the unit cell dimensions, the 
continuum description of graphene should breakdown. 

For this vortex solutions, the gauge field is linked with 
the angular momenta, relative to an axis perpendicular to 
the graphene sheet, of (p. Further, recall that the "mag- 
netic field" associated with this type of vortex solution 
vanishes and, therefore, the energy associated with the 
vortex also vanishes. 

For a vortex with a constant fo , equation (|52l) simpli- 
fies to 



1 



9v 



A 



0+ = O. (54) 



The solutions of this equation are the modified Bcssel 
functions / and K of argument Aqt for particular com- 
binations of the angular momenta associated with tp^ (/>+ 
and (j)- - see appendix lAl for details. 

The gauge model has, at least, one normalizable zero 
mode state of the Dirac equation. Therefore, the ?7&(1) 
and ?7j(l) (g) [/{,(!) gauge models can accommodate elec- 
tron fractionalization without the presence of external 
electromagnetic fields. More, given the vortex solution 
and relation ((53)) . besides fractionalization, the gauge 
models also incorporates flux quantization associated 
with the "magnetic field" . The phase-shifts coming from 
the Bohm-Aharonov effect associated with the vortex are 
connected with the component of the angular momenta 
of if along an axis perpendicular to the graphene layer. 



A. Zero modes and Conformal Invariance 
Breakdown 

The Sturm-Liouville form of the zero mode equation is 
found by setting = F/ y^, with z ~ Aq r, in equation 
([5^ . which becomes 



4 



F 



(55) 



where 



( m — I 
V J 



This equation is equiva- 



lent to a one-dimensional Schrodingcr eigenvalue prob- 
lem with a potential l/z^ which expresses invariance un- 
der scale change. The conformal symmetry is broken 



by the ultraviolet physics associated with defects/unity 
cell scales. Given that ly^ > 0, the potential strength 
for the Schrodinger problem is above the Breitenlohner- 
Freedman bound^ and the corresponding quantum me- 
chanical model is free from instabilities, i.e. the zero 
mode state does not collapse. In other contexts, for ex- 
ample in ultracold atomic physics, the violation of this 
bound gives rise to the Efimov effect^. 

We observe the connection between the fermionic 
Sturm-Liouville equation and the dynamics of fermion 
fields in a supergravity Anti-de-Sitter (AdS) background 
~ see, for example?^. In this description, the metric em- 
bodies conformal invariance leading to a potential 
associated, in our case, with the vortex solutions. The 
mass term of the fermionic field in the corresponding 
supergravity action contains the factor i^"^ (see e.gi^). 
Within this framework, the required short-range regu- 
larization (see the appendix) could be performed at the 
expense of introducing a dilaton field coupled to grav- 
ity, deforming the AdS metrio^i^. This suggests that 
the Maldacena conjecture of the AdS/CFT (conformal 
field theory) dualit y^^'^^ , may well provide fresh insights 
to graphene physics. For example, the vector mass gap 
could be a consequence of breaking exact symmetries in 
holographic 10-dimensional backgrounds, encoding mass 
gaps for the fermionic field and vector fiuctuations^. 



VI. RESULTS AND CONCLUSIONS 

In this paper we discuss gauge theories for graphene. 
The building of the gauge models starts assuming that 
graphene dynamics can be described by fermion fields 
together with a complex scalar field (p. The field p re- 
sumes the self-interaction of the carbon background and 
the mean fermionic self interaction. After exploring the 
global symmetry of the most general lagrangian, exclud- 
ing derivative-like couplings, the corresponding gauge 
models are investigated. 

The gauge models are compatible with a gap for 
fermion, vector and scalar fields. The mass gaps are 
generated via an Higgs mechanism. Further, the mass 
gaps associated with each gauge model are connected 
in different ways, which opens for the possibility of ex- 
perimentally distinguish between the models. Indeed, as 
claimed in section IIII Dl for example, changing the con- 
centration of impurities in graphene, one can change the 
fermionic gap and check how the scalar and vector gaps 
adjust and. in this way, check which of the gauge mod- 
els, if any, reproduce the graphene results. Furthermore, 
in what concerns the fermionic gap, within the models 
considered here, the mass gap is a dynamical quantity 
associated with the field (p. 

The gauge models have finite energy vortex solutions. 
Therefore, phenomena like the flux quantization of the 
"magnetic fleld" , in association with topological defects 
of the carbon structure, and/or Bohm-Aharonov type ef- 
fects become possible within the description of graphene 
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by gauge models. Several types of vortex solutions where 
discussed and, in general, the gauge field is linked with 
component of the angular momenta, along an axis per- 
pendicular to the graphene plan, of For this type 
of solution the phases associated with Bohm-Aharonov 
type effects are a measure of the Lp angular momenta. 

Finally, we have investigated the two dimensional 
Dirac equation for the vortex solutions. A generalization 
of the sublattice symmetry was discussed and we showed 
that all gauge models have normalizable zero modes to- 
gether with a non- vanishing fermionc gap. Fractional- 
ization is then possible in all the gauge models consid- 
ered here. Within this theoretical background one ex- 
pects that fractional quantum Hall effect can take place 
in graphene in connection with the zero mode solutions, 
even when there are no external electromagnetic fields. 
See also the discussion at the end of section H] 

The models investigated are potentially useful to de- 
scribe graphene. Indeed, they unify, under the same dy- 
namical principle, several features of the graphene and 
predict others. The models have multiple parameters 
whose values should be found by reproducing experimen- 
tally known graphene properties. We are currently en- 
gaged in performing such an investigation and will report 
the results elsewhere. 



Appendix A: Zero Modes of the Dirac Equation 

Let us discuss the solutions of equation ([54| . 



1 



b+ = Q. (Al) 



Introducing the adimensional distance z ~ Aq r, after 
multiplying this equation by one gets the following 
differential equation 



,2 



m —n I + 1 



= 



(A2) 



whose solutions are the modified Bessel functions I±,^{z) 
and K^{z), seo^ for definitions, where 



with its first derivative having a similar functional be- 
havior. Therefore, modulo the behavior for small r, in 
principle, setting (f>^ proportional to K^, the differential 
equation (j54p is solved and the spinor is normalizable. 
For small r 



K,{z) 



(A5) 



and the spinor diverge. 

The divergence of can be resolved as described in^ 
in connection with the potential — with /3 > 0. A 
short distance cut-off ro is introduced and the "poten- 
tial" is replaced by its value at tq. However, if (/>+ can 
be made regular at the origin, being given by (jSip 
will diverge near the origin. Due to this short distance 
divergence, is not normalizable unless the model has 
a minimal distance beyond which the continuum descrip- 
tion of graphene no longer makes sense. The dimensions 
of the unit cell provide such an infrared cut-off. Remem- 
ber that, in graphene, the carbon atoms are separated by 
a « 1.42 A and for r <a one can set 4>±{r) sa 4>±{a) and 
have a continuous and normalizable spinor. 

We call the reader attention that the introduction of 
a short distance, i.e. ultraviolet, cut-off for the fermions 
fields does not change the results of section IIV Al In- 
deed, the Dirac zero modes give no contribution to the 
equations of motion associated with the vortex solution 
of section IIV Al In this sense, the fiux quantization of 
the gauge field and the Bohm-Aharonov effect discussed 
there, i.e. the topological properties of the vortex solu- 
tion, are independent of the Dirac spinors. On the other 
hand, the requirement that the Dirac spinors are normal- 
izable is a necessary condition to have electron fraction- 
alization; see, for example, the discussion in the work of 
Hon et ati^. 

The singular behavior at r = seem to be an indi- 
cation of a "hole" at the center of the unit cell. This 
" hole" is a topological obstruction and is at the origin of 
the topological properties of the model analyzed in the 
present work. 
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m n 

9v . 



(A3) 



Note that v can be a real number. Of I±u{z) and K^{z) 
functions, only the last one tends to zero as r ^ +oo. 
Indeed, in this limit 



K,iz) 



(A4) 
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